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Abstract 

We present the applications of variational-wavelet ap- 
proach for the analytical/numerical treatment of the effects 
of insertion devices on beam dynamics. We investigate the 
dynamical models which have polynomial nonlinearities 
and variable coefficients. We construct the corresponding 
wavelet representation for wigglers and undulator magnets. 

1 INTRODUCTION 

In this paper we consider the applications of a new nume- 
rical-analytical technique which is based on the methods 
of local nonlinear harmonic analysis or wavelet analysis to 
the treatment of effects of insertion devices on beam dy- 
namics. Our approach in this paper is based on the general- 
ization of variational-wavelet approach from [l]-[8], which 
allows us to consider not only polynomial but rational type 
of nonlinearities [9] . We present solution via full multires- 
olution expansion in all time scales, which gives us expan- 
sion into a slow part and fast oscillating parts. So, we may 
move from coarse scales of resolution to the finest one for 
obtaining more detailed information about our dynamical 
process. In this way we give contribution to our full solu- 
tion from each scale of resolution or each time scale. The 
same is correct for the contribution to power spectral den- 
sity (energy spectrum): we can take into account contri- 
butions from each level/scale of resolution. Starting from 
formulation of initial dynamical problems (part 2) we con- 
struct in part 3 via multiresolution analysis explicit repre- 
sentation for all dynamical variables in the base of com- 
pactly supported wavelets. Then in part 4 we consider fur- 
ther extension of our previous results to the case of variable 
coefficients. 



2 EFFECTS OF INSERTION DEVICES 
ON BEAM DYNAMICS 

Assuming a sinusoidal field variation, we may consider ac- 
cording to [10] the analytical treatment of the effects of 
insertion devices on beam dynamics. One of the major 
detrimental aspects of the installation of insertion devices 
is the resulting reduction of dynamic aperture. Introduction 
of non-linearities leads to enhancement of the amplitude- 
dependent tune shifts and distortion of phase space. The 
nonlinear fields will produce significant effects at large be- 
tatron amplitudes such as excitation of n-order resonances. 
The components of the insertion device vector potential 
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used for the derivation of equations of motion are as fol- 
lows: 

Ax = cosh{kxx) cosh{kyy) sm{ks)/{kp) (1) 
Ay = k,j: smh(kxx) smh{kyy) sm(ks)/{kykp) 



with kl + ky 



k"^ — (27r/A)^, where A is the period length 
of the insertion device, p is the radius of the curvature in 
the field Bq. After a canonical transformation to betatron 
variables, the Hamiltonian is averaged over the period of 
the insertion device and hyperbolic functions are expanded 
to the fourth order in x and y (or arbitrary order). Then we 
have the following Hamiltonian: 
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[Pxiklx^ + ky) - 2klpyxy\ 



We have in this case also nonlinear (polynomial with de- 
gree 3) dynamical system with variable (periodic) coeffi- 
cients. After averaging the motion over a magnetic period 
we have the following related equations 
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3 WAVELET FRAMEWORK 

The first main part of our consideration is some varia- 
tional approach to this problem, which reduces initial prob- 
lem to the problem of solution of functional equations at 
the first stage and some algebraical problems at the sec- 
ond stage. Multiresolution expansion is the second main 
part of our construction. Because affine group of trans- 
lation and dilations is inside the approach, this method 
resembles the action of a microscope. We have contri- 
bution to final result from each scale of resolution from 
the whole infinite scale of increasing closed subspaces Vj'. 
...V^2 C K-1 C Vb C Vi C ^2 C .... The solution 
is parameterized by solutions of two reduced algebraical 
problems, one is nonlinear and the second are some linear 
problems, which are obtained by the method of Connec- 
tion Coefficients (CC)[11]. We use compactly supported 
wavelet basis. Let our wavelet expansion be 
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(4) 
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If Cjk = for j > J, then f{x) has an ahernative ex- 
pansion in terms of dilated scaling functions only f{x) = 
S cji(pji{x). This is a finite wavelet expansion, it can 

be written solely in terms of translated scaling functions. 
To solve our second associated linear problem we need 
to evaluate derivatives of f{x) in terms of (p{x). Let be 
= d'^(pi{x)/dx". We consider computation of the 
wavelet - Galerkin integrals. Let /"^{x) be d-derivative of 
function f{x), then we have /"^{x) — J^e^iVgi^)' 
values fi{x) can be expanded in terms of (p{x) 



{x) = ^A„(^m(x), (5) 

Xm = / ip'l{x)iprn{x)Ax , 



where A„i are wavelet-Galerkin integrals. The coefficients 
A„i are 2-term connection coefficients. In general we need 
to find {di > 0) 



A' 



did2---dj 



CO 



(6) 



For Riccati case we need to evaluate two and three connec- 
tion coefficients 

/oo 
^''^{x)^1-{x)dx, (7) 
-oo 

OO 



According to CC method [11] we use the next construc- 
tion. When N in scaling equation is a finite even positive 
integer the function ipix) has compact support contained in 

[0,7V-1]. For a fixed triple {di,d2,dj,) only some Kf^^'^'' 
are nonzero: 2-N<t<N~2, 2 - N < m < 
N-2, \i-m\ < N-2. There are Af = m^-9N + 7 
such pairs {£, m). Let A'^^'^^'^^ be an M-vector, whose com- 
ponents are numbers A' 



did2d3 



Then we have the first re- 
duced algebraical system : A satisfy the system of equa- 
tions (d = di + c?2 + ds) 

y^^did2d3 _ 2l-djydid2d3 ^g-j 
P 

By moment equations we have created a system of AI + d+ 
1 equations in M unknowns. It has rank M and we can ob- 
tain unique solution by combination of LU decomposition 
and QR algorithm. The second reduced algebraical system 
gives us the 2-term connection coefficients {d = di + d2)'- 

AKd.d2 ^2^~d^d,d2^ A, ,^ = ^apa,_2^+p (9) 

p 

For nonquadratic case we have analogously additional lin- 
ear problems for objects (6). Solving these Unear problems 



we obtain the coefficients of reduced nonlinear algebraical 
system and after that we obtain the coefficients of wavelet 
expansion (4). As a result we obtained the explicit time so- 
lution of our problem in the base of compactly supported 
wavelets. On Fig.l we present an example of basis wavelet 
function which satisfies some boundary conditions. In the 
following we consider extension of this approach to the 
case of arbitrary variable coefficients. 

4 VARIABLE COEFFICIENTS 

In the case when we have the situation when our problems 
(2), (3) are described by a system of nonlinear (rational) dif- 
ferential equations, we need to consider also the extension 
of our previous approach which can take into account any 
type of variable coefficients (periodic, regular or singular). 
We can produce such approach if we add in our construc- 
tion additional refinement equation, which encoded all in- 
formation about variable coefficients [12]. According to 
our variational approach we need to compute only addi- 
tional integrals of the form 



6„ (i)((^i)'^K2"i - fci)(<^2)''H2"t- fc2)dx, (10) 



where 6^ {t) are arbitrary functions of time and trial func- 
tions (/?!, (/?2 satisfy the refinement equations: 



[t) = ^ aikipi{2t - k) 



(11) 



If we consider all computations in the class of compactly 
supported wavelets then only a finite number of coefficients 
do not vanish. To approximate the non-constant coeffi- 
cients, we need choose a different refinable function (^3 
along with some local approximation scheme 

{Btf){x) FiAnM'^'t - k), (12) 

qGZ 

where F(^k are suitable functionals supported in a small 
neighborhood of 2^^k and then replace bij in (10) by 



Bgbijlt). In particular case one can take a characteristic 
function and can thus approximate non-smooth coefficients 
locally. To guarantee sufficient accuracy of the resulting 
approximation to ( |l0| ) it is important to have the flexibility 
of choosing (p^ different from ipi,ip2- In the case when D 
is some domain, we can write 



b^3{t) \d= J2 hj{t)XDi2't~k), 



(13) 



0<fc<2* 



where xd is characteristic function of D. So, if we take 
f4 ~ Xd, which is again a refinable function, then the 
problem of computation of ( [ic| ) is reduced to the problem 
of calculation of integral 

H{ki,k2,k3,k4) = H{k)= [ ^i{2H-ki)- 
ip3{2H - k2)(pf' (Tt - k3)ip'^^2H - ki)dx (14) 



The key point is that these integrals also satisfy some sort 
of refinement equation [12]: 

2-l^li/(fc) -^62fc-^ifW, fi = di+d2. (15) 

This equation can be interpreted as the problem of comput- 
ing an eigenvector. Thus, we reduced the problem of ex- 
tension of our method to the case of variable coefficients to 
the same standard algebraical problem as in the preceding 
sections. So, the general scheme is the same one and we 
have only one more additional linear algebraic problem by 
which we can parameterize the solutions of corresponding 
problem in the same way. 

On Fig. 2 we present approximated configuration and on 
Fig. 3 the corresponding multiresolution representation ac- 
cording to formula (4). 




Figure 1; Basis wavelet with fixed boundary conditions 
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Figure 2: Approximated configuration 
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Figure 3: Multiresolution representation 

ogy and Symplectic Scales', New Applications of Nonlin- 
ear and Chaotic Dynamics in Mechanics, 31,101 (Kluwer, 

1998) . 

[3] A.N. Fedorova and M.G. Zeitlin, 'Nonlinear Dynamics of 
Accelerator via Wavelet Approach' , CP405, 87 (American 
Institute of Physics, 1997). 
Los Alamos preprint, physics/9710035. 

[4] A.N. Fedorova, M.G. Zeitlin and Z. Parsa, 'Wavelet Ap- 
proach to Accelerator Problems', parts 1-3, Proc. PAC97 2, 
1502, 1505, 1508 (IEEE, 1998). 

[5] A.N. Fedorova, M.G. Zeitlin and Z. Parsa, Proc. EPAC98, 
930, 933 (Institute of Physics, 1998). 

[6] A.N. Fedorova, M.G. Zeitlin and Z. Parsa, Variational Ap- 
proach in Wavelet Framework to Polynomial Approxima- 
tions of Nonlinear Accelerator Problems. CP468, 48 (Amer- 
ican Institute of Physics, 1999). 
Los Alamos preprint, physics/990262 

[7] A.N. Fedorova, M.G. Zeitlin and Z. Parsa, Symmetry, 
Hamiltonian Problems and Wavelets in Accelerator Physics. 
CP468, 69 (American Institute of Physics, 1999). 
Los Alamos preprint, physics/990263 

[8] A.N. Fedorova and M.G. Zeitlin, Nonlinear Accelerator 
Problems via Wavelets, parts 1-8, Proc. PAC99, I6I4, 1617, 
1620, 2900, 2903, 2906, 2909, 2912 (IEEE/ APS, New York, 

1999) . 

Los Alamos preprints: |physics/9904039|, bhvsics/9904040l 
physics/9904041, physics/9904042, |physics/9904043| , pliy- 
sics/9904045, physics/9904046, physics/9904047. 

[9] A.N. Fedorova and M.G. Zeitlin, Los Alamos preprint 



3hysics/0003095 



[10] A. Ropert, CERN 98-04. 

[II] A. Latto, H.L. Resnikoff and E. Tenenbaum, Aware Techni- 
cal Report AD910708, 1991. 

[12] W. Dahmen, C. Micchelli, SIAM J. Numer. Anal, 30, 507 
(1993). 



